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Abstract — By introducing sliding mode predictive control (SMPC) 
techniques, a controller with passivity of a class of discrete-time 
uncertain singular systems is proposed in this paper. A new 
switching surface function is designed by taking the singular 
matrix into account, thus the resulting sliding mode dynamics is a 
full-order uncertain singular time-varying delay system. Due to 
feedback correction and receding horizon optimization, the 
influence of uncertainty can be compensated in time, strong 
robustness to matched or unmatched uncertainties is possessed. In 
addition, chattering of sliding mode control can be eliminated by 
predictive control method. Simulation result is given to illustrate 
the validity of the proposed approach. 
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I. 



INTRODUCTION 



Singular systems, also referred to as descriptor systems, 
generalized state-space systems or differential-algebraic 
systems, provide convenient and natural representations in the 
description of economic systems, power systems and circuits 
systems. So a great number of problems about singular systems 
were studied (see [l]-[4]). Time delays are frequently the main 
causes of instability and poor performance of systems, and are 
encountered in a variety of engineering systems such as 
chemical processes, nuclear reactors, biological systems and so 
on. 

Sliding mode control (SMC) is a popular control approach 
for systems containing uncertainties or unknown disturbances, 
as the controllers can be designed to compensate for the 
uncertainties or disturbances. For a broad class of systems, this 
kind of control is particularly appealing due to its ability to 
deal with non-linearities, time-variance, as well as uncertainties 
and disturbances, in a direct manner in the face of modeling 
imprecisions (see [5]-[7]). The first step in SMC is to define a 
sliding surface. The second step is to design the control law in 
such a way that any state outside the sliding surface is driven to 
reach the surface in finite time and stay here. 

However, it is well known that chattering is a flaw for SMC, 
to reduce chattering usually results in decreasing robustness of 
the closed-loop systems. So a SMPC strategy is proposed to 
overcome this problem in this paper. Model predictive control 
(MPC) is a control technique which permits to cope with a 
constrained system providing an optimal control strategy. The 
concepts of prediction and receding horizon that generalized 
prediction control provides, can improve the performance in a 
reaching mode, and this is one of the goals that the SMPC 
strategy achieve (see [8]-[15]). Another one is the capability of 
controlling processes with large time delays and high 



controllability ratio. Moreover, the problem of the 
implementation of a SMC when the state is not accessible can 
be solved with the predictive strategy. 

There are many works based on SMPC control techniques. 
For example, [8] applied SMPC in a solar air conditioning 
plant, SMPC for nonlinear systems based on lazy learning is 
investigated in [9], robust model predictive control with 
integral sliding mode in continuous-time sampled-data 
nonlinear systems is studied in [10], [11] investigated model- 
based predictive networked control systems. But they are 
continuous case. Due to the widespread using of digital 

Implementations, discrete-time systems are very common 
in real plants compared with continuous-time ones. Besides, 
stable continuous-time systems may become unstable after 
being discretized. Therefore, it is necessary to design 
controllers for discrete-time systems directly. [12] investigated 
a SMPC algorithm for a class of discrete-time n -joint rigid 
robotic manipulator systems. Stability analysis for a triangular 
discrete-time nonlinear system with SMPC approach studied in 
[13]. 

In this work, we will address the problems of passivity- 
based SMPC for a class of uncertain discrete-time singular 
time-varying delay systems. We will pay particular attention to 

the singular matrix -k in the design of a new switching surface 
function, which leads to a full-order uncertain singular time- 
delay system for describing the sliding mode dynamics. Then 
construct a newly Lyapunov function to derive a delay- 
dependent sufficient condition in the form of linear matrix 
inequality (LMI), which guarantees that the sliding mode 
dynamics is robustly passive. A numerical example illustrates 
the validity of the proposed approach. 

Notation: Throughout this paper, " denotes the n - 

dimensional Euclidean space, R is the set of all mxn real 

matrices; V - ) denotes P as a positive definite 

(positive semi-definite) symmetric matrix. » x \\ refers to the 

Euclidean norm of the vector x ; •* denotes the identity matrix 
with proper dimension. Matrices, if not explicitly stated, are 
assumed to have compatible dimensions. The symmetric terms 
in a symmetric matrix are denoted by * . 



II. 



PROBLEM FORMULATION 



Consider the following time-varying delay discrete-time 
singular system: 
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Ex(k +1) = Ax(k) + A d x(k-d(k)) + B ro co(k) 

+ Bu(k), (1) 

z(k) = Cx(k) + C d x(k-d(k)) + D w a>(k), 
x(k) = <?>(k),ke[-d,Q,] 

, x(k) e R" . t , . . u(k) € i? m • t . t , • t 

where v ' is the state, v ' is the control input, 

z(k)eR p . ., . , . . w(k)ei? q . ., 

v ' is the control output, v ' is the 

disturbances, *■ ' represents the time-varying delay and 

satisfied — — *■ ' ~ , ™\ ) i s the initial function; 

e K may be singular and it is assume that 

rank(E) = r<n ; A,A d ,B,B m ,C,C d md D m ^ ^ 

matrices with appropriate dimensions. 

We recall the following definitions for the discrete-time 
singular time-delay system. 

JEx(k+l) = Ax(k) + A d x(k-d(k)), 
[x(k) = <f>(k),te[-d,0.] 

Definition 1. The system (1) is said to be : 

(i) regular if ^ 'is not identically zero, 

(ii) causal if it is regular and 
deg(det(zE - A)) = rank(£) 

(iii) stable if it is regular and "*■ ' ^ ' , where 



+d(k)T] T (k) 



W 1 



R 1 ^ Y 2 W]Tj(k), 



(4) 



and 



(2) 



D(0,1) 



expresses the interior of the identity circle which 



center is the origin, 

p(zE-A) = {z\d t(zE -eA) = 0} . 

Definition 2. The uncertain singular time-delay system (1) 
with *• ' is said to be robustly passive if there exists a 

scalar ' such that 



2^ r (k)z(k)>- r ^ r (k)#) 



s=0 



(3) 



Lemma 3. (Integral Inequality) For any matrices > ' 
v 2 ' and a time-varying delay *■ > , then 

- X Z T (s)E T REm 



s=k-d(k) 



<^ T (/c) 



E^+Y^E E T Y 2 +Y 1 T E E T W 

* -E T Y 2 +Y 2 T E -E T W 

* * n 



Tj(k) 



where <?(s) = x(s+l)-x(s) 

?j(k)=[f(k) f(k-d(k))0) T (k)] 

III. MAIN RULES 

In this section, the problems of regularity, causally, stability 
and passivity for singular system (1) are investigated. 
Furthermore, a controller is given by SMPC techniques which 
make the closed-loop system is robust stability and passivity. 

A. Passivity Analysis 

The following theorem shows that the system (1) 

( U W — U) | s passivity can be guaranteed if there exist some 
matrices satisfying a certain LMI. 

Theorem 1. Given a scalar ' , the discrete-time singular 

system (1) is regular, causal, stable and robustly passive, if 

p,Q ™h z 



there exist positive matrices 
d> o> i anc j 2 satisfying 



and 



matrices 



11 * 12 


1 13 


^ 


^15 


A T P 


w 

J. 22 


1 23 


dY 2 T 


dAjZ 


AjP 


« * 


1 33 


dW T 


dBlZ 


BlP 


« * 


* 


-dZ 








« * 


* 


* 


-dZ 





« * 


* 


* 


* 


-P 



<o, 



(5) 



where "12 ~~ " Q- 



av§ 



V F 11 = A T i?S r + SR T A-E T PE + Y*E + E T Y l 

+ (d 12 +l)Q, 
Y 12 = A T i?Sj + SR T A d + E T Y 2 -F/£, 
^i3 = SR r B CJ + A r i?S:-C r , 
V 15 =d n (A-E) T Z, 

Y 22 = AjflSj + S d i? T A d - Y 2 T E - E T Y 2 - Q, 
^¥ 23 = S d R T B a + A^RSl-C T d , 
y - - y i + BlBSl + SlR T B n - D n - Dl , 

where R is any full column matrix that satisfied ^ "■ ~ . 
Proof: Choose a Lyapunov function candidate as 
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v(fc) = £v,(fc) 



where 



V 1 (/c) = x T (k)E T PEx(k), 
V 2 (k)= £ x T (s)Qx(s), 



s=/c-d(k) 
-d+1 fc-1 



v 3 (fc)= E E * t ooq*(s)> 



0=-d+2s=k-l+0 



V 4 (fc)= Z E z T (s)E T ZEz(s). 



e=-d+is=k-i+e 



By Lemma 3, we have 

AV^k) = x T (k + l)E T PEx(k + l)-x T (k)E T PEx(k), 
AV 2 (k)<x T (k)Qx(k)-x T (k-d(k))Qx(k-d(k)) 

+ § * T (s)Q*(s), 

s=k+l-d 

AV 3 (/c) = (d-d)x r (/c)Qx(/c)- j x r (s)Qx(s), 

s=k-d+l 

AV 4 (/c) = dz r (/c)E r ZEz(/c)- E z T (s)E T ZEz(s), 



s=k-d 



Then 



AV(k)<?i T (k){ 



A T 
gr 



(A-£) T 



P[A 4 B a ] 



dR[(A-E) A d Bj 



+ 



-E J PE+Q 

* -Q 

* * o 

E T Y 1 + Y*E -YfE + E T Y 2 E T W 
* -E T Y 2 -Y 2 T E -E T W 



W T 



dZ- 1 ^ Y 2 W]}Tj(k), 



where 



Tj(k)=[ X (k) X T (k-d(k)) 0) T (k)] T , 

z(s) = x(s + l)-x(s). 

Notice that ^ "■ ~ " , we have 

2x T (k +1)E T R(S T x(k) + S T d x(k-d(k)) 

+S>(fc)) = 0. 

Therefore, 

AV(Jc) - 2e/ (k)z(k) - r« r (k)o(fc) 
<<f(k)Q<f(k), 

where 



n = 



fu 


1 12 


1 13 


< 


^15 


A r P 


* 


a 22 


- 1 23 


d7 2 r 


<Z 


AjP 


* 


* 


- 1 33 


dW T 


dBlZ 


B^P 


* 


* 


* 


-dZ 








* 


* 


* 


* 


-dZ 





* 


* 


* 
1 13 


* 

y 13 + c 


* 


-P 






- 1 23 


^ 23 +C 


d> 








- 1 33 


B^I 


+ S** 1 


B„. 



(6) 



(7) 



(8) 



(9) 



Sum both sides of (8) from to -K gives rise to 

fc-1 k-l 

V(JC)-V(0 }2^o) T (k)z(k)-r^a) T (k)o)(k)<0 (10) 



Under the zero initial condition, we have ^ ' and 

V(K) > ^ so ^ guarantees (3) 

Next we prove the singular system (1) is regular and causal. 
Since ran K(E,)-r<n ^ mere must exist two i nve rtible 
G,HeR" 



matrices 



such that 
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GEH 


= 











Thenj? = G r 


"0" 

o 


, GAH = 


Ai A2 

A 21 A 22 


? 


T PG T = 




p 12 

p 

L x 22_ 


,G~ T Y 1 H = 


Y Y 

± n 12 

y y 

/21 ^22_ 



s(/c) = GEx(k) - Y,G(A-E)x(s) 



(13) 



S 12 
c- ^n <0 A (d-d+l)Q>0 

Since lx and v — y ^ , we get 

:= A T RS T + SR T A-E T PE + Y^E + E J Y 1 < 0. 

Pre- and post-multiplying " ^ " by " and " , respectively. 
Then 



Al 2 <^S T 21 + S 21 cD r A 22 < 0, 



(11) 



i.e., 



A2 



A 

^ 2 is nonsingular. Otherwise, supposing J ' 22 is singular, 

there must be a non-zero vector ' which ensure that 

A C — 

^ 2 ' . And then we can conclude that : 

t T (Al 2 ®S T 21 + S 21 ® T A 22 K = 0, 
A 

and this contradicts (11). So ^ 2 is nonsingular. Then the 
singular system (1) is regular and causal. 

Finally, we consider the stability of nominal case of system 

(1) with ^ ' . Following the similar procedure as used 
above, we obtain 



AV(k)<ff T (k) x ¥ff(k), 



(12) 



where #)4/(fc)/(fc-d(fc))f and 



where A and^> are system matrices defined in (1), the matrix 
" is to be chosen such that * J ' " is nonsingular. 

Remark 1. Notice that the discrete-time sliding surface 
function designed in (13) is different from that previous 

literatures, since the singular matrix E is considered in this 

paper. This enables to avoid some difficulties caused by E in 
deriving the sliding mode dynamics subsequently. 

C. Design of Sliding Mode Prediction Model 

In order to implement an MPC, a model of the plant is used 
to predict the future plant outputs. This prediction is based on 
past and current values of the input and the output of the plant. 

The sliding mode prediction model is constructed 

s m (k +1) = GAx(k) + GA d x(k -d(k)) + GBu(k) 

-^G(A-E)x(s) + Zs(k), (14) 

s=0 

where £ is a designable parameter, which satisfies < E, < 1 . 

D. Design of Control Law 

In practice, because of time-variance, non-linearity or 
disturbances, SMPC will inevitably exist errors, therefore, the 
model output will not the same as real output. A desirable way 
to solve such problem is feedback correction 

h(*+l) = s m {k+l) + o(s{k)-s m {k)), (15) 

where °" is a weight coefficient, i.e., weighted feedback 
correction. The effect of feedback correction will reduce with 
the decreasing of °" . When suitable °" is chosen, feedback 
correction can compensate the model and improve control 
performance. From the viewpoint of practice, the appropriate 

range of <? is 0<cr<1 . 



V F = 



Fii 


T 12 


< 


^15 


A T P 


* 


± 22 


dY 2 T 


dAjZ 


AjP 


* 


* 


-dZ 








* 


* 


* 


-dZ 





* 


* 


* 


* 


-P 



s(k) = s(k)-s m (k) 



Let 
Then (15) reduces to 



for the sake of clarification. 



Noxwy xr) ax (12) i\xnXisa * < " , xrjsv 

Ay(/c)<-l min (-^)|^(/c)| 2 <0. 

Thus, the stability of the nominal system (1) is proved. The 
following, a SMPC controller will be designed. 

B. Sliding Mode Function 

Defining sliding mode function as 



S in (k+l) = S m (k+l) + CTS(fc) 

= GAx(k) + GA d x(k-d(k)) + GBu(k) 

k 

-J]G(A-E)x(s) + £s(k) + crs(k) 

s=0 

= h(k) + GBu(k), 

where 

h(k) = GAx(k) + GA d x{k - d(k)) 

k 

~Y J G{A-E)x{s) + £s(k) + as(k) . 

s=0 



(16) 
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Now, performance index is given as 

J=(s m (/c+l)-s r ) 2 +Au 2 (/c), 
s 

where r is a sliding mode reference value. 



fe-i 



(17) 



Since the control objective is to keep states on the sliding 
surface, the desired sliding mode reference value should be 

s r =0 

Therefore, the performance index (17) can be reduced to 

J=s 2 m (k + 1) + Au 2 (k), (18) 

where A is a weight coefficient, which adjusts optimized 

index ^ and the control signal. 

According to equation (16), performance index (18) can be 
rewritten as 

J = (h(k) + GBu(k)) 2 +Au 2 (k). (19) 

The solution of minimizing (19) gives control signal *■ ' . 
By setting the partial derivative of ^ to zero, i.e., 
OJ I Oil — U ^ solving the resulting equation, the optimal 
solution of *■ ' is 



u{k) = . 



h(k)GB 
U + {GBf 
u(k) 



(20) 



Now, the control input v ' for the closed-loop system (1) is 
obtained. 

E Robust Stability Analysis 

From performance index (19), one can see that *• ' 

affects -' less with the decreasing of ^ . For the sake of 
simplifying analysis, according to optimal control theory (see 

[13]), it is sound to suppose "- = " , i.e., the case when control 
signal is not optimal. Therefore, control law (20) reduces to 

u(k) = -h(k)(GBy\ (21) 

Consider the singular system (1) and sliding mode function 
(13), it follows 



+GBu(/c-l)-^G(A-£)x(s) 

s=0 

+£s(/c-l) 
= GEx(k)-GB a o)(k-l) 

k-l 

-jG(A-£)x(s) + ^-l) (23) 

s=0 

Accordingly, 

s(k) = s(k)-s m (k) 

= GB a <o(k-l)-£s(k-l). (24) 

Therefore, (22) turns to 

s(k+l) = -4s(k)-&s(k) + GB a a>(k) 
= -£s(k)-oiGB a a>(k-l) 
-ts(k-l)\ + GB m a>(k) 
= ^[-s(k) + o-s(k-l)] 

+GB a [a>(k)-oe>(k-l)] 
= {M+GB a N, (25) 

where M =-s(k) + crs(k-l),N = co(k)-<j(a(k-l). 

Theorem 2. If the change rate of disturbances is bounded, 
i.e., the following inequality holds 



\o}(k)-CTO}(k-i)\< 



M 



GB, 



(26) 



s(/c + l) = GEx(/c+l)-^G(A-£)x(s) 

s=0 

= GAx{k) + GA d x{k-d{k))-h(k) 
+GB a cD(k)-^G(A-E)x(s) 



= GB a co(k)-Zs(k)-cjs(k). 
According to (14), 

s m (k) = GAx(k-l) + GA d x(k-d(k)-l) 



(22) 



where " is a positive constant, then the closed-loop system (1) 
which is constructed by (13) and (20) is regular, causal, 
robustly stable and passive. 

Proof: Consider the characteristic polynomial of -<" t 

-1+ crz' 1 = 0. (27) 

Obviously, the root of equation (27) is z ~ & • Because 
0<ct<1 ; M is stable. 

While ' , ' is stable. Namely, ' ° 

such that ' ' ' when ° . 

According to (27), I I - ^ , then 

\s(k+l)\ = \^M + GB a N\ 

<\EM\ + \GB\\N\< s + u. 

I II (9|| • 

Consequently, the practical sliding mode motion of the 
closed-loop system will converge to a £ vicinity of sliding 
surface and stay on it subsequently. In addition, because the 
stability of sliding surface has been guaranteed by (13), the 
closed-loop system (1) is robustly stable with the control (20), 
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the regularity, causally and passivity can be get from Theorem 
1 directly. 

Remark 2. It is well known that chattering is a flaw for 
SMC, to reduce chattering usually at the cost of decreasing 
robustness of the closed-loop systems. In SMPC, because of 
feedback correction and receding horizon optimization, the 
influence of uncertainties may be discovered in time, the 
control signal can be adjusted immediately to prevent system 
states cross sliding surface, hence chattering will be avoid. 
Thus the combination of SMC and predictive control can do 
great help to controller design. 

IV. SIMALATION RESULTS 

In this section, an example is given to illustrate the 
advantage of proposed method. 

Example 1. Consider the uncertain singular time-varying 
delay system (1) with the following parameters: 



-0.3 -0.2 0.3 
0.5 0.3 0.5 
0.9 0.4 0.9 



optimization. The usefulness of the proposed method has been 
verified by the numerical results. 



"1 0" 




1 


,A= 








0.1 
0.8 
0.5 
0.5 

0.5 







0.2 



" 




"0.5" 


0.3 


,B = 


0.5 







0.5 



C = [0.1 0.2 0.1], 



"0.2" 




"0.2" 


0.1 


A,= 


-0.2 


0.3 




0.4 



Set Lr — r> A -with ■& — 1 for computational simplicity, 

the initial state of system is ^ ' 6 5 - 4 Lr= 0-02, ct = .95. 
From (2 1), we get the controller. The simulation results of the 
closed-loop system (1) with the proposed method in this paper 
are illustrated in Figs. 1. 

It is well known that chattering is a flaw for SMC. In 
SMPC, because of feedback correction and receding horizon 
optimization, the influence of uncertainties may be discovered 
in time, the control signal can be adjusted immediately to 
prevent system states cross sliding surface, hence chattering 
phenomena is eliminated with SMPC method and the close- 
loop system (1) is stable. 

V. CONCLUSION 

This paper has investigated the problems of robust passivity 
analysis and SMPC of uncertain singular time-varying delay 
systems. The major theoretical findings are as follows. First, 
the delay-dependent sufficient condition in the form of LMI 
has been established so as to ensure that the sliding mode 
dynamics is stable and robustly passive. Then the desired 
control law is constructed by using special sliding mode 
prediction model, feedback correction and receding horizon 



J 










i ' 




XjM 

x 2 (k) _ 

x 3 (k) 


V\ 
















i \ 


/ 
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